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Abstract 

Minimax solutions are weak solutions to Cauchy problems involv- 
ing Hamilton-Jacobi equations, constructed from generating families 
quadratic at infinity of their geometric solutions. We give a complete 
description of minimax solutions and we classify their generic singu- 
larities of codimension not greater than 2. 
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1 Introduction 

Hamilton-Jacobi equations play an important role in many fields of math- 
ematics and physics, as for instance calculus of variations, optimal control 
theory, differential games, continuum mechanics and optics. 

Let us consider a Cauchy Problem involving a Hamilton-Jacobi equa- 
tion. For small enough time t the solution u is classically determined using 
the characteristic method. Although u is initially smooth, there exists in gen- 
eral a critical time beyond which characteristics cross. After this time, the 
solution u is multivalued and singularities appear. Therefore, the problem 
how to extract "true solutions" from multivalued solutions naturally arises. 



"To appear in Journal of Geometry and Physics. 
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In 1991 Marc Chaperon proposed in [3] a geometric method to construct 
weak solutions to Hamilton- Jacobi equations, called minimax solutions (see 
also [13]). Their definition is based on generating families quadratic at in- 
finity of Lagrangian submanifolds. The minimax solutions have the same 
analytic properties as the viscosity solutions, namely existence and unique- 
ness theorems hold. However, minimax solutions are in general different from 
viscosity solutions. 

This paper is organized as follows. The first part is a survey on minimax 
solutions. In the second part we study the small codimension generic sin- 
gularities of minimax solutions. Namely, we prove that the singularities of 
codimension not greater than 2 of minimax solutions and viscosity solutions 
are the same. This result may be useful to find examples in which minimax 
solutions have physical meaning. 

The survey part of the paper is mainly based on a talk I have given for the 
"Trimester on Dynamical and Control Systems" at SISSA-ICTP in Trieste. 

1 thank Andrei Agrachev for its kind invitation. 

2 Minimax of functions quadratic at infinity 

Let E be the space R n and / : E — > R a smooth function having a finite 
number if critical points. We suppose / quadratic at infinity, that is / is a 
non-degenerate quadratic form outside a compact set. For A G R, we set 

E x := G E : /(£) < A} . 

Let us fix a G R, big enough for that (—a, a) contains all the critical values 
of /. Then we set E ±ao := E ±a . 

Let us denote by Q\ the restriction mapping E — > E x , inducing the homo- 
morphism : H*(E, E~°°) — > H*(E X , E~°°) between the relative reduced Z- 
homology groups. The quotient space E/E~°° is homeomorphic to § 9 , where 
q denotes the index of the quadratic part of / at infinity. Hence, the only 
non-trivial reduced homology group is H q (E, E~°°) ~ Z. Fix a generator 7 
of this group. 

Definition. The minimax of / is the real number 

minmax(/) := inf{A G R : ^7 7^ 0} . 
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The infimum defining the minimax exists and it is finite, indeed g^y = 
7 7^ and Q*_xl = whenever A > a; it does not depend on the choice of the 
generator 7. Moreover, the minimax of a function is a critical value, since by 
definition the topology of the sublevel sets E x changes when A crosses the 
minimax value. 

Since the minimax level is stable under small deformation of the func- 
tion (see j2]), we may assume without loss of generality that / is generic. 
Genericity conditions mean that all the critical points of / are of Morse type 
and all its critical values are different. Therefore, every critical point has an 
index, which is the number of negative signs in the Morse normal form of the 
function at the critical point. 

Morse Theory (see 0) describes how the topology of the sublevels sets 
E x changes when A changes. Namely, E x is diffeomorphic to E** , provided 
that [A, fj] does not contain critical values. On the other hand, if [A, fj] con- 
tains only one critical value, then E^ retracts on the space obtained from E x 
attaching to its boundary a cell of dimension equal to the index of the critical 
point realizing the critical value. Hence, every generic function on E defines 
a cell decomposition of the space and there is a 1 : 1 correspondence between 
its cells and the critical points of the function. 

Let us recall that the incidence coefficient of two cells in this decompo- 
sition is the degree of the restriction of the attaching map from the higher 
dimension cell to the lower dimension cell (see [B]). The incidence coefficient 
[£ : rf[ of two critical points £ and 77 of / is the incidence coefficient of the 
corresponding cells. Note that [£ : 77] 7^ implies that ind(£) = ind(ry) + 1 
and /(£) > f(rj). For example, the two critical points of /(£) = £ 3 — e£ have 
non-zero incidence coefficient for e > 0. 

The incidence coefficient of pairs of critical points leads to a natural par- 
tition of the critical set £ of / in the following way. Denote by C\ a pair of 
critical points realizing the minimum of the set 

- m ■ t,v e £, [£ :v] ¥> 0} . 

Then define by induction Cj + i from C\,...,Ci as a pair of critical points 
realizing the minimum of the set 

{/(0 - f(v) ■ t,V G £ \ (Ci U • • • U d), [£ : rj] ? 0} . 

In this way, we decompose the critical set into the disjoint union of pairs Ci 
and a set F which does not contains incident critical points, i. e. [£ : rj] — 
for every £,77 G F. 
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Definition. Two critical points are coupled if they belong to the same pair 
Cj in the preceding decomposition; a critical point is free if it belongs to F. 

In [2 J we proved the following result. 

Theorem 1. Every generic function quadratic at infinity has exactly one 
free critical point, and its value is the minimax. 

3 Geometric and multivalued solutions 

In this section we introduce the symplectic framework for the geometric and 
multivalued solutions of Hamilton-Jacobi equations. Let X be a closed man- 
ifold of dimension n, n : T*X — > X its cotangent bundle, endowed with the 
standard symplectic form dX, where A is the Liouville's 1-form. In local co- 
ordinates T*X = {x,y}, we have Tc(x,y) = x and A = y dx. A dimension 
n submanifold of T*X is said to be Lagrangian whenever dX vanishes on 
it. Two Lagrangian submanifolds are Hamiltonian isotopic if there exists 
a flow generated by a Hamiltonian field, transforming one into the other. 
Hamiltonian isotopies transform Lagrangian submanifolds into Lagrangian 
submanifolds. A generating family of a Lagrangian submanifold L C T*X is 
a smooth function S : X x M. k —>■ R such that 



Given a generating family S, the following operations give rise to new 
generating families T of the same Lagrangian submanifold: 

1. Addition of a constant: T(x; £) = S(x; £) + C for C G R; 

2. Stabilization: T(x; £, 77) = S(x; £) + Q(r]), where Q is a non-degenerate 
quadratic form; 

3. Diffeomorphism: T(x;rj) = S(x; £(:c, 77)), where (x; 77) 1— >■ (x, 77)) is 
a fibered diffeomorphism with respect to the coordinate x. 

Two generating families are equivalent if one can obtain one from the 
other by a finite sequence of the preceding operations. A generating family 
5* is quadratic at infinity (gfqi) if there exists a non-degenerate quadratic 
form Q such that S(x;^) = <3(£), whenever |£| is big enough (uniformly in 




d e s=o = max. 
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Theorem 2 (Sikorav— Viterbo, [10J — [12j ) . Every Lagrangian subman- 
ifold ofT*X, Hamiltonian isotopic to the zero section {(x;0) : x G X}, 
admits a unique gfqi modulo the preceding equivalence relation. 

Remark. The Theorem still holds in the case of non-compact manifolds, 
provided that the projection of the Lagrangian submanifold into the base is 
1 : 1 outside a compact set. 

A Lagrangian submanifold L C T*X is exact if the Liouville's 1-form is 
exact on it. In this case, we can lift L to a Legendrian submanifold in the 
contact space J X X of 1-jets over X. This lifting projects into a wave front 
in the space J°X ~ X x R = {x, z} of 0-jets over X, defined up to shifts in 
the z-direction. If 5* is a gfqi of L, then the corresponding wave front is the 
graph of S 



Graphs of equivalent gfqi are equal up to a translation along the 2-axis. The 
points for which the natural projection "forgetting z" of the graph into X 
is not a fibration form a stratified hypersurface in X, provided that L is 
generic. 

We may now define generalized solutions to Hamilton- Jacobi equations. 
Let Q be a closed manifold. We consider the following Cauchy Problem on 
Q involving a Hamilton-Jacobi equation: 



The Hamiltonian H : R + x T*Q ->• R is supposed of class <?f 2 on R+ x T*Q 
and continuous at the boundary; the initial condition w : Q — > R is assumed 
to be of class c £ 1 . 

Let M := Q x R be the space-time manifold and T*M = {t,q;r,p} its 
cotangent bundle (endowed with the symplectic form dp Adq + dr A dt) . We 
consider the flow $ : R + x T*M — > T*M, generated by the Hamiltonian 

: T*M R defined by Jf := r + H, and the submanifold 



a := { ((0, q; -H(0, q, du (q)), du (q)) : q G Q} C J^~\0) C T*M . 



{{x,S{x]S)):dzS{x;O = 0} . 




Definition. The geometric solution of {CP) is the submanifold 




5 



It turns out that the geometric solution to (CP) is an exact Lagrangian 
submanifold, contained into the hypersurface J^ 7_1 (0) and Hamiltonian iso- 
topic (for finite times) to the zero section of T*M (see [2]). Hence, the 
Sikorav-Viterbo Theorem guarantees that it has a unique gfqi modulo the 
equivalence relation. Up to a suitable constant, we may assume that the 
graph of the gfqi restrained at t = is equal to a. We denote by S such a 
gfqi. 

Let us consider now the projection pr : T*M — > T*Q, defined locally 
by pr(t, q;r,p) := (q,p). We call geometric solution at time t of (CP) the 
exact Lagrangian submanifold L t := pro$*(a) C T*Q, which is an isochrone 
section of the global geometric solution. It is easy to check that L t is Hamil- 
tonian isotopic to the zero section, so, by the Uniqueness Theorem, its gfqi 
is S t (q;0 ■■= S(t,q;£). 

Definition. The graph of S (resp., St) will be called the multivalued solution 
(resp., at time t) to {CP). 

Remark. It is possible to construct global generating families of geometric 
solutions as follows. The action functional J pdq — Hdt is a global generating 
family, whose parameters belong to an infinite dimension space. By a fix 
point method, proposed by Amann, Conley and Zehnder, one can obtain a 
true generating family (with finite dimensional parameters), see 

4 Minimax and viscosity solutions 

In this section we introduce minimax solutions to Hamilton-Jacobi equations, 
and we discuss their properties and relations with viscosity solutions. Let S 
be the gfqi of the geometric solution to {CP). For every fixed point (t, q) in 
the space time, the function S t>q is quadratic at infinity, so we may consider 
its minimax critical values, defined in section El 

Definition. The minimax solution to (CP) is the function defined by 

u(t,q) : = minmax(St )g ) . 

Remark. For a Cauchy Problem (CP) on a non-compact manifold, as for 
instance for Q = K n , we define the minimax solution as follows. Let (H n ) n be 
a sequence of smooth Hamiltonians with compact support in the variables q. 
Then, at any Cauchy Problem defined by H n and the given initial condition, 
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we may associate its minimax solution U n , denned as above. By definition, 
the minimax solution of (CP) is the limit of the solutions U n for n — > oo. 

Theorem 3 (Chaperon, |4j). The minimax solutions are weak solutions 1 
to (CP), Lipschitz on finite times, which not depend on the choice of the 

gfqi- 

We end this section discussing the relations between minimax and vis- 
cosity solutions. Viscosity solutions to Hamilton-Jacobi equations have been 
introduced by Crandall and Lions in [H], following the generalization proposed 
by Kruzhkov of the Hopf 's formulas. 

Theorem 4 (Joukovskaia, |8j). The minimax and the viscosity solutions 
of (CP) are equal, provided that the Hamiltonian defining the equation is 
convex or concave in the variables p. 

Izumiya and Kossioris constructed in [7] the viscosity solution beyond its 
first critical time. Actually, they proved that the Lagrangian graph of the 
viscosity solution is not always contained into the geometric solution of the 
equation. Hence, minimax and viscosity solutions are in general different. 

5 Characterization of minimax solutions 

In order to describe a useful characterization of minimax solutions (first 
presented in [2J), we start introducing some standard notations concerning 
unidimensional wave fronts. 

Let J°R ~ R 2 = {(q, z)} be the space of 0-jets over R and vr : J°R -> R 
the natural fibration (q, z) i— > q. A wave front in J°R is the projection of 
an embedded Legendrian curve in the contact space J 1 !^ ~ R 3 = {(q,z,p)} 
under the projection 7Ti : (q,z,p) h- > (q,z). The only singularities of generic 
wave fronts are semicubic cusps and transversal self-intersections. Two em- 
bedded Legendrian curves in J^-R are Legendrian isotopic if there exists a 
smooth path joining them in the space of the embedded Legendrian curves. 
In this case their projections are also called isotopic. A front is long if it is 
the graph of a smooth function outside a compact set of R; it is flat if its 

1 A weak solution to (CP) is a continuous and almost everywhere derivable function, 
which solves in these points the Hamilton-Jacobi equation; moreover, its restriction at 
t = is equal to the initial condition. 
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tangent lines are nowhere vertical. A section of a flat front is a connected 
maximal subset which is the graph of a piecewise c <9 1 function; a branch is 
a smooth section. The generic perestroikas occurring to flat wave fronts un- 
der isotopies among flat fronts are the following, illustrated in figure ^ cusp 
birth/ death, triple intersection, cusp crossing. 



Figure 1: Generic singularities of flat fronts under Legendrian isotopies. 

We come back now to the Cauchy Problem (CP) in the case Q — R. 
We assume that its geometric solution is a generic Lagrangian curve. A 
multivalued solution £ of (CP) at a generic time is a flat wave front, isotopic 
to {(q,0) : q G Up to replace our Cauchy Problem with a sequence of 
problems approximating it, we may assume also that £ is a long wave front, 
and that it is the graph of a gfqi, denoted by S, of the geometric solution. 
Hence, every branch of the multivalued solution is realized by a critical point 
of S tt g (depending on q as a parameter). 

Definition. The index of a branch is by definition the index of its corre- 
sponding critical point minus the index of the quadratic part of 5*. 

The index of a branch does not depend on the choice of the gfqi. The non- 
compact branch of £ has index 0, as well as the minimax section. Moreover, 
the index changes by +1 (resp., —1) passing through a positive cusp 2 (resp., 
negative cusp). A double point of £ is called homogeneous whenever it is the 
intersection of two branches having the same index. 

2 A cusp is positive (resp., negative) if, following the orientation of the front, we pass from 
a branch to the other according to the coorientation defined by the positive z direction. 




3/2 



3/2 
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The partition of the critical set of a generic function quadratic at infinity 
we give in section |21 leads, by Theorem El to a decomposition of the multi- 
valued solution into the graph fi of the minimax solution and the pairs Xi 
of coupled sections, formed by the sections corresponding to coupled critical 
points of the gfqi: 

E = //u|Jjq . 

i 

Indeed, when q runs on E, every pair of coupled critical points of St, q moves 
on the front along two sections, forming a closed curve. 

The curves X, are homeomorphic to a disc boundary, they have exactly 
two cusps and no self-intersections. They are almost everywhere smooth and 
only continuous at homogeneous double points of the multivalued solution 
(see 0). 

Any homogeneous double point of an oriented multivalued solution defines 
in it a closed connected curve. Such a curve is called a triangle whenever it 
has exactly two cusps. The double point is the vertex of the triangle. Let T 
be a triangle of S. Fix an arbitrary small ball centered at the vertex of T, 
intersecting only the two branches containing the double point. We denote 
by E — T a long flat front equal to E \ T outside the ball and smooth inside 
it (see figure |2J). 




Figure 2: The front E — T. 

Definition. A triangle T of E is vanishing if there exists a Legendrian isotopy 
joining E and E — T among long flat fronts without homogeneous triple 
intersections. 

The following Theorem (see [2]) provides an effective method to simplify 
recursively a given multivalued solution. After a finite number of steps, we 
get a front, which is actually the graph of the minimax solution to the initial 
problem. 
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Theorem 5. Let £ be a multivalued solution and fiUXiU- ■ -UXjy its decom- 
position into minimax section and coupled sections. Then £ is smooth or has 
a vanishing triangle T among the curves Xi. In this case, the decomposition 
of £ — T is induced by that of £. In particular, the minimax of E and E — T 
are equal outside the ball centered at the vertex of the vanishing triangle. 

Example. Consider the multivalued solution E depicted in figure El The 
homogeneous double points B, C, D, F, G and / are not verteces of triangles. 
The triangles T% and Tl of verteces E and L are vanishing, while those 
of verteces A, H and M are not. Indeed, to cut off Ta and T M we must 
pass through a self-tangency, while to cut off T H we must pass through a 
homogeneous triple intersection at F. The minimax of E is equal, outside two 
arbitrary small balls centered at E and L, to the minimax of £' := Tj—Te—Tl. 
Now, B and / define two vanishing triangles Tb and TV; £' —Tb—Tj is smooth, 
so it is the minimax section of E outside two arbitrary small balls centered 
at B and /. 




Figure 3: The multivalued solutions E and £'. 
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6 Singularities of minimax solutions 



In this section we classify the generic singularities of small codimension aris- 
ing in minimax solutions to Hamilton- Jacobi equations. This classification 
is made with respect to the Left-Right equivalence, fibered on the time di- 
rection, defined below. Let us consider two functions /, g : K x M n -> 1. We 
denote by t the coordinate in R and by q = (qi, . . . , q n ) the coordinate in R n . 

Definition. Let a, b be two points of the space-time R x R n . The germ of 
/ at a and the germ of g at b are said to be equivalent if there exists two 
diffeomorphism germs (p and ip such that ip o / o ip^ 1 = g and if is fibered 
with respect to the time axis: (p(t,q) = (T '(t) , Q (t , q)j. 

A singularity of a function germ for this equivalence relation is its equiva- 
lence class. A Cauchy Problem on Q = R n is said to be generic if its geometric 
solution is generic as Lagrangian submanifold; its minimax solution is also 
called generic. 

In jH] , Joukovskaia proved that the singular set of any generic minimax so- 
lution (formed by the points where the solution is not is a closed stratified 
hypersurface, diffeomorphic at any point to a semi-algebraic hypersurface. 

Theorem 6. Let u be a generic minimax solution and (t, q) a point belonging 
to a stratum of codimension c = 1 or 2 of its singular set. Then the germ of 
u at (t, q) is equivalent to one of the map germs in the table below. 



c 


Normal form 


1 


\<li\ 




2 


min{ \qi 


It} 


2 


min{y 4 - tY 2 + 


qi Y :Y G R} 



The images of these normal forms (in the case n = 2) are depicted in 
figure |U Note that these singularities are stable. These map germs are also 
the normal forms of the generic singularities (of codimension not greater than 
2) of viscosity solutions (whose classification has been done by Bogaevski in 

my 

Joukovskaia studied in jS] the generic singularities of minimax functions, 
defined as minimax levels of Lagrangian submanifold's gfqi which are not 
necessarily geometric solutions of Hamilton-Jacobi equations. To prove The- 
orem |U1 we will show that some singularities in Joukovskaia's list do not arise 
in minimax solutions. In order to do this, let us recall her Theorem. 
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Figure 4: Small codimension generic singularities of minimax solutions. 

Theorem 7 (Joukovskaia, |8j). Any generic minimax function on R x R" 
is equivalent, at any point in a codimension 2 stratum of its singular set, 
to one of the map germs listed in Theorem (for c = 2) or to one of the 
following map germs: 

(a) (t; 5i, . . . ,q n ) i-> max {t, 

(b) (t; q u . . . , q n ) ^ min {{q^, max{-|gi|, t}}. 

The images of the map germs (a) and (b) are shown in figure El 




Figure 5: Further generic singularities of minimax functions. 
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Proof of Theorem® We may assume n = 1 without loss of generality. Let 
us consider a multivalued solution having a singularity (a) or (b). In both 
cases, there are three 0-index branches crossing at the singular point (settled 
at the origin). 

Fix an arbitrary small neighborhood U of the origin in K; then for every 
small enough time t < 0, the isochrone multivalued solution at this time 
takes, over the corresponding section of U, the configuration illustrated in 
the left part of figure 13 moreover, all the other sections of the front's decom- 
position are smooth over U. The minimax section \i has three jumps in A, 
B and C, shrinking to the homogeneous triple point at the origin as t goes 
to 0. 

Actually, the jumps A and C belong to the same section, denoted by 
a in the figure. By Theorem El we may recursively eliminate the vanishing 




Figure 6: Multivalued solution configurations near the singular point. 



triangles in the multivalued solution changing neither the minimax section 
nor the other pairs of coupled critical points. Hence, after a finite number of 
such operations, the front is composed by the minimax solutions, the section 
a and its coupled section a*, see the right part of figure |H1 

This front is not isotopic to the trivial front {(q, 0) : q 6 R}. Hence, mul- 
tivalued solutions have neither singular points of type (a) nor (b). Theorem 
IHlis now proven. □ 

Remark. In other words, both singularities (a) and (b) imply the existence of 
a cycle in the graph associated to a decomposition of an isochrone multivalued 
solution of the Cauchy Problem. This is impossible, since all the graphs 
associated to admissible decompositions of isochrone multivalued solutions 
are trees, due to Theorem El 
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